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Abstract 



We study the complexity of computing average quantities related to spin systems, such 
as the mean magnetization and susceptibility in the ferromagnetic Ising model, and the average 
dimer count (or average size of a matching) in the monomer-dimer model. By establishing 
connections between the complexity of computing these averages and the location of the 
complex zeros of the partition function, we show that these averages are #P-hard to compute. 
In case of the Ising model, our approach requires us to prove an extension of the famous 
Lee-Yang Theorem from the 1950s. 

1 Introduction 

1.1 Background 

Many natural computational problems in combinatorics, statistics and statistical physics can be 
cast in the following framework. We are given as input a graph G = (V, E) which implicitly 
defines a set £1 = Q(G) of combinatorial structures, or configurations (such as matchings in G, 
or /c-colorings of its vertices). A weight function w : £1 — > R + assigns a positive weight to every 
element a 6 fi, giving rise to a probability distribution n(a) = w(a)/Z; here the normalizing 
factor Z := Y2aen w(a) is called the partition function, and is typically not known. (While it is 
easy to compute the weight w(a) for any given a, computing the sum Z is usually hard. Note 
that the size of is typically exponential in the size of G.) 

The final ingredient is a non-negative function, or observable, f : Q — > R + U {0}, so that f(a) 
is also easily computable for each a G $7. Our goal is to compute the average of / with respect 
to 7r, i.e., 



As illustration we present two classical examples, both of which we will develop extensively in 
our later results. 

Example 1 : The Ferromagnetic Ising Model. Here the configurations are assignments of spin values 
{+, -} to the vertices of G, i.e., £1 = {+, — } v . The weight of a configuration a is 



where d(a) is the number of disagreements in a (i.e., the number of edges {u, v} e E with 
a(u) / cr(v)), and p(a) is the number of vertices v £ V with a(v) = +. The model has two 
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parameters: the edge potential (3, satisfying < (3 < 1, which governs the strength of the inter- 
action between neighboring spins; and the vertex activity A > 0, which specifies the tendency 
for spins to be +. The probability distribution ir(a) = wi(a)/Zj is the familiar Gibbs distribu- 
tion, and Zi := Zi (G, (3, A) is the associated partition function. Note that when f3 < 1, this 
distribution favors agreement between neighboring spins. Similarly the distribution favors '+' 
(respectively '— ') spins when A > 1 (respectively when A < 1). 

An important observable here is the magnetization p(a), which is just the number of +-spins 
in a. Its average, the mean magnetization, is a fundamental quantity in statistical physics: 



Zi 

Other widely studied averages include the mean energy (d) (the average size of the cut between 
+-spins and — spins) and the susceptibility x '■= (p 2 ) - (p) 2 (the variance of the magnetization). 

□ 

Example 2: Matchings, or the Monomer-Dimer Model. The configurations Q are all matchings 
(independent sets of edges) in G. The weight of a matching a is 

w M (a):=X u ^H 7e , (2) 

where u(a) is the number of unmatched vertices (monomers) in a. The parameter A > is 
the vertex weight (or monomer activity), while for each edge e G E, j e is an edge weight (or 
dimer activity). The Gibbs distribution it (a) = wm(&) / Zm is a natural weighted distribution 
on matchings, and the partition function Zm '■= Zm(G, {"f e }eeE, A) is the weighted matching 
polynomial of G. 

A natural observable here is u(a), the number of unmatched vertices (or monomers). Note 
that (|V| - (u))/2 is just the average size of a (weighted) matching in G (or equivalently, the 
average number of dimers). □ 

Observe that the Ising model partition function Zi may be written as a polynomial in A (actually 
a bivariate polynomial in A and (3): 

\v\ 

Zi = ^a k X k , where a k = £ f3 d ^\ 

k=0 a:p(a)=k 

The mean magnetization then becomes 

w _£^*_^, (3) 

where V denotes the differential operator \ Similarly, the mean energy and susceptibility \ 
can be written 

(d> = /3^-; x-={p)-{p) = -2ir~{~zT) ' ( ) 

For matchings, the partition function Zm is nothing other than the matching polynomial 

\v\ 



Zm = ^2 "fcA 
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k=0 



2 



where a& = Yla-u(a)=k Tleecr is a weighted sum over matchings with k unmatched vertices. 
The average number of monomers is then 

VZm 

(«> = — , (5) 

where V again denotes the differential operator A^. 

Equations (3)-(5), which express averages as the ratio of some derivative of the partition 
function to the partition function itself, are in fact no accident; they are a consequence of the 
fact that the Gibbs distribution takes the form w[a) = exp(—H(a))/Z, where the Hamiltonian 
H(a) is a sum of natural observables. 

The subject of this paper is the computational complexity of computing natural averages 
such as (3)-(5). While the complexity of computing partition functions has been widely studied 
in the framework of Valiant's class #P of counting problems (see, e.g., [9-12, 17, 19]), we are not 
aware of any corresponding results for the exact computation of averages. In the approximate 
setting, by contrast, it is well known that (at least for the wide class of self-reducible problems, 
which includes all the examples above) approximate computation of the partition function is 
polynomial time equivalent to sampling (approximately) from the Gibbs distribution ir [24]; 
and sampling from ir clearly allows us to approximate averages to any desired accuracy (this 
is because the observables one is concerned with in these situations admit a priori absolute 
bounds which are polynomial in the input size). 

What if we are interested in exact computation? It is tempting to argue that computing an 
average as in, say, (3) is at least as hard as computing the partition function Zj, because (3) is 
a rational function and thus by evaluating it at a small number of points we could recover the 
numerator and denominator polynomials by rational interpolation. Since the partition function 
is #P-hard in almost all cases of interest (including Zj and Zm above at all but trivial values of 
the parameters), we would be done. 

The problem with this argument is that, viewed as polynomials in the variable A, Zi and its 
derivative T>Zi may have common factors (equivalently viewed as polynomials in the complex 
variable A, they might have common zeros); and in this case we are clearly not able to recover Zi 
by rational interpolation. Indeed it seems hard a priori to rule out the possibility that non-trivial 
interactions between Zi and its derivative could conspire to make the average much easier to 
compute than Zi itself. Thus we are naturally led to the following question: 

Question: Is it possible for the partition function Z and its derivative to have common zeros? 1 

If the answer is no, then we will be able to conclude that computing the average is as hard as 
computing Z itself, and thus #P-hard in all interesting cases. This is the main goal of this paper. 



1.2 Contributions 

The question of common zeros actually turns out to be a deeper issue of wider interest in statis- 
tical physics and complex analysis. The study of the zeros of the partition function dates back 
to the work of Lee and Yang in 1952: the famous Lee-Yang Circle Theorem [26] proves the remark- 
able fact that the zeros of the ferromagnetic Ising partition function Zj always lie on the unit 
circle in the complex plane. This classical theorem, which has since been re-proved many times 
in different ways [2, 28, 35], was developed initially as an approach to studying phase transi- 
tions, but has since spawned a more global theory connected with the Laguerre-Polya-Schur 
theory of linear operators preserving stability of polynomials. (See the Related Work section 
below.) 

1 Note that a common zero of Z and VZ corresponds to a repeated zero of Z, so this question is equivalent to the 
question of whether Z may have repeated zeros. 
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Somewhat surprisingly, despite much activity in this area, the question of the location of 
the zeros of the derivative VZj (or equivalently of repeated zeros of Zj itself) remains open. 
The main technical contribution of this paper is to resolve this question as follows. 

Theorem 1.1. Let G = (V, E) be a connected graph, and sup-pose < j3 < 1. Then the zeros of the 
polynomial VZi(G, (3, A) (in X) satisfy |A| < 1. 

Since the Lee-Yang Theorem says that all zeros of Zj satisfy |A| = 1, Theorem 1.1 immedi- 
ately implies that Zj and VZj have no common zeros. 

Remark. The restriction that G be connected is needed: there exist disconnected graphs for 
which the conclusion of the theorem does not hold. A simple example is a graph consisting of 
two isomorphic disconnected subgraphs. For the same reason we require (3 < I. We also note 
that standard facts from complex analysis (in particular, the Gauss-Lucas theorem) imply that 
the zeros of VZj lie in the convex hull of those of Zi, and hence within the closed unit circle. 
The content of Theorem 1.1 is that they must lie in the interior of the circle. This refinement is 
of course crucial for our application. 

Before moving on, let us briefly mention our approach to the proof. We actually prove a 
more general result concerning the zeros of the multivariate partition function Zj(G, (3, { X v } v& y )) 
see Theorem 3.2 in Section 3. (The Lee-Yang Theorem itself is also often stated in multivariate 
form.) Our proof is based on a delightful combinatorial proof of the Lee-Yang Theorem due 
to Asano [2], which begins with the empty graph (which trivially satisfies the theorem) and 
builds the desired graph G by repeatedly adding edges one at a time; by a careful induction 
one can show that the Lee-Yang property is preserved under each edge addition. Our proof 
follows a similar induction, but the argument is more delicate because we are working with the 
more complicated polynomial VZj rather than Zj. In particular, in the inductive step we need 
to evoke a non-trivial correlation inequality due to Newman [28] . 

Our first computational complexity result follows as an almost immediate corollary of The- 
orem 1.1. 

Theorem 1.2. For any fixed < (3 < 1 and any fixed A / 1, the problem of computing the mean 
magnetization of the Ising model on connected graphs is #P-hard. Moreover, the problem remains #P- 
hard even when the input is restricted to graphs of maximum degree at most IS., for any fixed A > 4. 

Note that in the case A = 1 the mean magnetization is trivially \ V\/2 by symmetry. Theo- 
rem 1.2 confirms that in all non-trivial cases, the problem of computing the fundamental aver- 
age quantity associated with the Ising model is as hard as it could possibly be. Furthermore, 
the result also holds for bounded degree graphs, which are relevant in the statistical physics 
setting. The result can also be extended to arbitrary ferromagnetic two-spin systems and to 
planar graphs: the details can be found in Appendix C. 

We also prove a similar (but slightly weaker) result for the susceptibility of the Ising model. 

Theorem 1.3. For any fixed < j3 < 1, the problem of computing the susceptibility of the Ising model 
on connected graphs, when A is specified in unary, is #P-hard. Moreover, the problem remains #P-hard 
even when the input is restricted to graphs of maximum degree at most A for any A > 3. 

Remark. The requirement that A be part of the input seems to be an artifact of the rational 
interpolation operations we use in our proof. In particular, our proof of Theorem 1.2 shows 
hardness for fixed A by "simulating" different values of A by suitably modifying the graph. To 
adapt this reduction approach to prove hardness for susceptibility (at fixed values of A) seems 
to require the polynomial time computation of magnetization as a subroutine. However, we 
conjecture that computing the susceptibility should be hard even for fixed values of A (including 
A = l). 
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We then proceed beyond the Ising model, and ask about the hardness of computing aver- 
ages in the monomer-dimer model (i.e., weighted matchings). A classical result of Heilmann 
and Lieb [23] establishes an analog of the Lee-Yang Theorem for the zeros of the monomer- 
dimer partition function Zm', however, Heilmann and Lieb also present examples of (connected) 
graphs G for which Zm has repeated zeros, so we cannot hope to prove an analog of Theorem 1 . 1 
in this case. On the other hand, Heilmann and Lieb show that if G contains a Hamiltonian path 
then all the zeros of Zm are simple. We are able to capitalize on this fact by adapting exist- 
ing #P-hardness reductions for Zm in such a way that the instances of Zm that appear in the 
reduction always contain a Hamiltonian path. Specifically, we present a reduction from the 
problem Monotone 2-SAT of counting satisfying assignments of a monotone 2-CNF formula 
to computing Zm in Hamiltonian graphs G. The reduction is an elaboration of Valiant's original 
#P-completeness proof for the permanent [36]. 

This leads to our third computational complexity result. 

Theorem 1.4. For any fixed A > 0, the problem of computing the average number of dimers (equiva- 
lently, average size of a matching) in the monomer-dimer model on connected graphs with edge weights 
in the set {1, 2, 3} is #P-hard. Moreover, the problem remains #P-hard even when the input is restricted 
to graphs of maximum degree at most A, for any A > 5. 

Remark. Note that our hardness result requires a small finite number (three) of different values 
for the edge weights. However, this requirement can be removed if G is allowed to have parallel 
edges; the theorem then holds for any single fixed non-zero edge weight (including the uniform 
case in which all edge weights are 1). 

1.3 Related Work 

The study of the location of zeros of the partition function was initiated by Yang and Lee [39] in 
connection with the analysis of phase transitions. In the follow-up paper [26], they instantiated 
this approach for the ferromagnetic Ising model by proving the celebrated Lee-Yang theorem 
on the location of zeros of the partition function and using it to conclude that the ferromagnetic 
Ising model can have at most one phase transition. The Lee-Yang approach has since become 
a cornerstone of the study of phase transitions, and has been used extensively in the statisti- 
cal physics literature: see, e.g., [2, 23, 28, 35] for specific examples, and Ruelle's book [32] for 
background. In a slightly different line of work, Biskup et al [4,5] studied a novel approach to 
Lee-Yang theorems for a general class of spin systems on lattice graphs using asymptotic expan- 
sions of the partition function. Zeros of partition functions have also been studied in a purely 
combinatorial setting without reference to the physical interpretation: see, for example, Choe 
et al [14] for a collection of such results about zeros of a general class of partition functions. 
Another important example is the work of Chudnovsky and Seymour [15], who show that the 
zeros of the independence polynomial of claw-free graphs lie on the real line. There have also 
been attempts to relate the Lee-Yang program to the Riemann hypothesis [29]. 

Lee-Yang theorems have also been studied in mathematics in connection with the theory of 
stability preserving operators. The main problem underlying this area is the characterization of 
linear operators that preserve the class of polynomials, called Q-stable polynomials, whose zeros 
lie in some fixed closed set 11 This research area has its origins in the work of Laguerre [25] and 
of Polya and Schur [31], and also has connections to control theory [16] and to electrical circuit 
theory [8]. It has also seen considerable recent activity, especially through the breakthrough 
results of Borcea and Branden, who completely characterize stability preserving operators for 
multivariate polynomials in various important settings [6,7]. Although the study of stability 
preserving operators is closely related to our problem, there is a crucial difference in that we 
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require our linear differential operator to not only preserve the stability of the partition function, 
but in fact to improve it, by restricting the possible locus of the zeros of the derivative to the open 
interior of the locus of the zeros of the partition function itself. 

In the statistical physics literature, we are aware of only two works which consider the multi- 
plicity of the zeros of the Ising partition function: Heilmann and Lieb [23] and Biskup et al [4,5]. 
In [23], a theorem similar to our Theorem 1.1 is proven in the special case when the underly- 
ing graph G has a Hamiltonian path and (3 is close enough to 1 (depending upon the graph 
G). Similarly, in the special case of the Ising model, the results of [4] imply our result but only 
when P is close to 0, and only in the special case of lattice graphs [3,5]. Note that neither of 
these results appears to be sufficient for the purposes of our hardness result. 

The classification of counting problems associated with partition functions (via so-called 
dichotomy theorems) has also recently been a very active area of research. For several interesting 
general classes of partition functions, these theorems characterize the partition function as be- 
ing either computable in polynomial time or #P-hard [9-12, 17, 19]. However, there appear to 
be no analogous results on the complexity of averages such as the magnetization. 

A related area that we do not deal with in this work is the problem of approximate count- 
ing. Recent progress in this area has shown that the complexity of approximating the partition 
function, as well as that of the related problem of approximate sampling, is closely related to 
the phase transition phenomenon [33,34,38]. However, it is not clear whether hardness results 
analogous to [33, 34] can be proven for the approximate computation of the magnetization. 

2 Preliminaries 
2.1 The models 

Let G = (V, E) be an undirected graph. The two models we will be concerned with are the 
ferromagnetic Ising model and the monomer dimer model, both of which have already been defined 
in Section 1.1. 

Ferromagnetic Ising model. Recall that in the ferromagnetic Ising model, a configuration a : V — > 
{+, — } is an assignment of +/— spins to the vertices of G. The model is characterized by an 
edge potential < f3 < 1, and a vertex activity A > 0. The weight function wi(a) defined in 
(1) induces a probability distribution over configurations with an associated partition function 
Zi(G, (3, A) := wi(a). We shall be concerned with the mean magnetization M(G, (3, A) := (p), 
which is the average number of +-spins in a configuration, and the susceptibility \ '■= {p 2 } — {p) 2 r 
which is the variance of the same quantity. As in (3)-(4), these quantities can be written in terms 
of the derivatives of Zi with respect to A. 

For our discussion of the zeros of Zj(G, (3, A), we will also need a generalization of the Ising 
model in which the vertex activities can vary across vertices of G. Suppose that the vertex ac- 
tivity at vertex v is z v . The weight of a configuration a is then defined as 



and the partition function is given by Zi(G, j3, (z v ) veV ) = s ^ JU wj{a). Consider the linear dif- 
ferential operator T>q defined as follows: 



Wl (a) := f3 d ^ ]"[ 



z. 



'V i 
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As in (3)-(4), we can then write the magnetization Mj(G, /3, (z v ) veV ) as 



M(G, ft, (z v ) veV ) 



Zr(G,f3,(z v ) veV ) 



Monomer-dimer model. Recall that in the monomer-dimer model, the configurations are match- 
ings of G. The model is characterized by edge weights 7 e > for every edge e in E and a vertex 
activity A > 0. The weight wm{p) of a matching a is as described in (2), and the associated 
partition function is defined by Zm(G, (7 e ) eg £ , A) := J2a w m{&)- 

The average number of monomers U(G, (7 e ) eg £ > ^) := ( u ) can ^ e written (as in (5)) in terms 
of the derivative of Z: in particular U (G, (7 e ) eg £ > ^) = V Zm{g{^)^ E \) ■ ^ ne avera S e dime* count 
D(G, ( , je) e& E ' ^) (equivalently, the average size of a matching) can be obtained from U by the 
simple relation 

g(g.(*W) = "-" ( °^ W) . 
where n is the number of vertices in G. 

Remark. In our definitions above, vertex activities are restricted to be positive real numbers. 
Although this is the physically (and computationally) relevant setting, in our proofs and in our 
discussion of Lee-Yang theorems we will need to work with vertex activities that are arbitrary 
complex numbers. The expressions for the quantities defined above still remain valid. 



2.2 Zeros of partition functions 

We first consider the location of the complex zeros of the partition function of the ferromagnetic 
Ising model. In a seminal paper Lee and Yang proved the following striking theorem [26]. 

Theorem 2.1 ([26]). Let G be any undirected graph and suppose < @ < 1. Then the complex zeros of 
Zm(G, j3, z), considered as a polynomial in z, satisfy \z\ = 1. 

Actually, Lee and Yang proved the following multivariate version of their theorem, the proof 
of which was later considerably simplified by Asano [2]. 

Theorem 2.2 ([2,26]). Let G = (V,E) be a connected undirected graph, and suppose < j3 < 1. 
Suppose (z v ) V £v is a set of complex valued vertex activities such that \z v \ > 1 for all v e V, and 
\z u \ > I for at least one u G V. Then Zj(G, f3, (z v ) veV ) ^ 0. 

Theorem 2.2 is readily seen to imply Theorem 2.1 by setting z v = z for all v 6 V. We now 
consider the partition function of the monomer-dimer model. In [22], Heilmann and Lieb stated 
the following result (see [23] for the complete proof). 

Theorem 2.3 ([22,23]). Let G = (V, E) be any graph, and ("y e )e&E be a collection of positive real edge 
weights. The complex zeros of Zm(G, (j e ) e ^E, z), considered as a polynomial in z, satisfy !ft (z) = 0. 
Further, if G contains a Hamiltonian path, all the zeros are simple. 

In [23], Heilmann and Lieb also gave examples of connected graphs G in which Zm has 
repeated zeros. This is in contrast to the Ising model, where, as we prove in Section 3, connect- 
edness is sufficient to ensure that the zeros are simple. 
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2.3 Rational interpolation 



In our hardness reductions, we will need a few well known facts about interpolation of rational 
functions. While it is clear that it is not in general possible to determine all coefficients of a 
rational function given its values at any number of points, this can be done if we impose a few 
simple conditions, as stated in the following theorem. 

Theorem 2.4 ([27]). Suppose R(x) = |£4 where gcd(p(x),q(x)) = 1 and both p(x) and q(x) are of 
degree n. Suppose p(x) and q(x) are polynomials of degree at most n satisfying 

for 2n + 2 distinct values xi, x%, . . . X2n+2- Then there is a constant c such that p(x) = cp(x) and 
q{x) = cq(x). 

Notice that given the evaluations at the points Xi one can write down a system of 2n + 2 
homogeneous linear equations for the 2n + 2 unknown coefficients of p and q. The theorem 
then guarantees that this system has rank exactly 2n + 1. Thus, since Gaussian elimination can 
be implemented to run in strongly polynomial time (see, e.g., [18]), a polynomial time algorithm 
for evaluating R immediately yields a polynomial time algorithm for determining some p and 
q satisfying the conditions of the above theorem. If we know at least one non-zero coefficient 
of p or q, we can then determine the proportionality constant c, and hence p and q also, in time 
polynomial in n. 



3 An extended Lee-Yang theorem 

In this section we prove Theorem 1.1, our extension of the classical Lee-Yang theorem. Let 
G = (V, E) be a connected graph with |V| = n and \E\ = m, with vertex activity z.- t at the ith 
vertex. When clear from the context , we will write Z{G) and M(G) for the partition function 
Zj(G, (3, (z v ) v& v) and the mean magnetization M(G, j3, (z v ) v& v) of the Ising model on G. In 
terms of the linear operator Vq defined in Section 2.3, we then have M(G) = T>gZ(G)/Z(G). 

For convenience, we will use the shorthand Y' = VqY (when G is clear from the context) 
in this section. Notice that this is slightly non-standard, as this shorthand is usually used for 
the actual derivative. In particular, when all the Z{ are equal to z, we have Y' = z^ with our 
notation. Also, observe that the operator Vq obeys the usual product rule: (Y1Y2)' = + 

In our proof, we will also need the following generalization of the partition function. We 
call an assignment of positive integer valued weights w : V — > Z + to the vertices of G legal if 
w(v) is at least equal to the degree of v, for all v G V. 

Definition 3.1. Let w be a legal collection of weights for G. The weighted partition function Z w (G) 
is then defined as 

z w (g) -.= ^ n < 6 ) 

<?£{+-} V v.a(v)=+ 

where, as before, d(a) is number of disagreeing edges in the configuration a. 

Notice that the multivariate Lee-Yang theorem (Theorem 2.2) holds also for the weighted par- 
tition function, since all the weights are positive integers and we are effectively just changing 
variables from z v to z™ . 

We will also need the following consequence of a correlation inequality of Newman [28], 
whose proof can be found in Appendix A. 
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Theorem 3.1 ([28, Theorem 3.2]). Let G be any graph, and let wbea legal collection of weights for G. 
Suppose < p < 1, and \z v \ > lfor allv^V are such that Z W (G) / 0. Then 2 



Here, 5ft (z) denotes the real part of z. 

In the special case of real valued activities, the above theorem is equivalent to the well known 
Griffiths inequality [21], which states the intuitive fact that in a ferromagnetic Ising model 
where all activities favor the + spin, the magnetization must be at least n/2. 

For ease of reference in the inductive proof, we give a name to the property we want to 
establish. Recall that when all the vertex activities are equal to z, the classical Gauss-Lucas 
theorem, together with the Lee-Yang theorem, implies that the zeros of the derivative T>g{Z{G) ) 
lie on or inside the unit circle. Our goal is to establish that they actually lie inside the unit circle. 
Accordingly, we use the following terminology: 

Definition 3.2 (Strict Gauss-Lucas property). A graph G = (V, E) has the strict Gauss-Lucas 
property (SGLP) if for every set of activities such that \z v \ > 1 for all v G V, and every < /3 < 1, 
one has T>gZ(G) / 0. The graph has the weighted strict Gauss-Lucas property (WSGLP) if for all 
legal weights w, V G Z W (G) / necessarily holds under the same conditions. 

Note that WSGLP easily implies SGLP: we simply choose w(v) = A for all v, where A is the 
maximum degree of G. From WSGLP, we then have that whenever \z v \ > 1 and < (3 < 1, 
V G Z W (G) = AV G Z(G) / 0, and hence V G Z{G) / 0. Thus Theorem 1.1 is implied by the 
following more general statement. 

Theorem 3.2. Every connected graph has the weighted strict Gauss-Lucas property, and hence also the 
strict Gauss-Lucas property. 

We now proceed to prove Theorem 3.2, using induction on the number of edges in the graph 
G. We first consider the base case of a connected graph with a single edge. 

Lemma 3.3 (Base Case). Let G be the graph consisting of a single edge connecting two vertices. Then 
G has the weighted strict Gauss-Lucas property. 

Proof. In this case we have Z W (G) = z^z™ 2 + f3(z^ + z™ 2 ) + 1 and therefore V G Z W (G) = 
(w\ + w 2 )z™ 1 Z2 2 + ^(wiz™ 1 + W2Z2 2 ), with wi,w 2 > 1. When \zi\ , \z 2 \ > 1, the latter vanishes 
only if 



For the inductive case, we require two operations: adding a new vertex to the graph, and 
merging two existing vertices. These operations are formalized in the following lemmas. 

Lemma 3.4 (Adding a vertex). Suppose G = (V, E) is a connected graph satisfying the weighted 
Gauss-Lucas property. Let ube a vertex not in V. Then, the graph G\ obtained by attaching the new 
vertex u to any vertex (say v±) of G also has the weighted Gauss-Lucas property. 

2 Recall that we are using here the slightly non-standard notation Z' w (G) = VgZ w (G), as described at the begin- 
ning of this section. 





which cannot hold since < f3 < 1. 



□ 
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Lemma 3.5 (Merging vertices). Suppose G = (V, E) is a connected graph satisfying the weighted 
strict Gauss-Lucas property. Consider any two vertices, say v\ and V2, in G that are not connected by 
an edge. The graph G\ obtained by merging v\ and V2 into a single vertex v (while making all the edges 
incident on v\ and v% incident on v) also has the weighted strict Gauss-Lucas property. 

Before proceeding with the proofs of the above lemmas, we show how to use them to prove 
Theorem 3.2. 

Proof of Theorem 3.2. We will prove by induction on m that any connected graph with at most 
m edges satisfies WSGLP. By Lemma 3.3, this statement is true when m = 1. Now suppose that 
the statement is true when m = k, and consider any connected graph G with k + 1 edges. 

In case G has a cycle, there exist vertices u and v such that the edge {u, v} can be removed 
from G to obtain a connected graph H. Since H has at most k edges, H satisfies WSGLP by 
the inductive hypothesis. Let v\ be a vertex not in G. By Lemma 3.4, the graph H U {{u, v\}} 
satisfies WSGLP. We can now merge v\ and v to obtain G, which therefore satisfies WSGLP by 
Lemma 3.5. 

In case G is a tree, there exists an edge {u, v} such that v is of degree 1. Again, we ob- 
tain a connected graph H with at most k edges by removing the edge {u, v}. By the inductive 
hypothesis, H satisfies WSGLP, and hence by Lemma 3.4, G does too. This completes the in- 
duction. □ 

Remark. Note that the proof of Theorem 3.2 given above holds also when the graph G is al- 
lowed to have parallel edges and self-loops. This will be useful in our extension to general 
two-state ferromagnetic spin systems in Appendix C. 

We turn now to the proofs of Lemmas 3.4 and 3.5, for which we will need the following 
additional lemma. 

Lemma 3.6. Let G be a connected graph. Fix any set S of vertices of G, and let Z+(S) denote the 
partition function restricted to configurations on the subgraph G — S, with all the vertices in S fixed to 
have spin +. Consider any set of vertex activities satisfying \z v \ > lforv EG — S. Then, for < (3 < 1 
and any set of permissible weights on the vertices of G,wehave Z+(S) / Oand?R. (Z+(S)' / Z+(S)) > 0. 
In particular, for any positive real a, we have Z+(S)' + aZ+(S) / 0. 

Proof. Observe that Z+(S) is proportional to the product of weighted partition functions on 
connected components of the graph G — S, where the activities on the vertices connected to S in 
these components (of which there is at least one in each component) have increased in magnitude 
by a factor of at least > 1. We can therefore conclude using Theorem 2.2 that Z+(S) / 
0. The second condition $l(Z+(S)' /Z+(S)) > then follows from Theorem 3.1 applied to 



We first prove Lemma 3.5, since its proof is somewhat simpler. 

Proof of Lemma 3.5. Consider any legal weight assignment on G\. If the weight of v in G\ is w v , 
we can write w v = w\ + u>2 such that the weight assignment giving weights w\ and W2 to v\ 
and v 2 respectively is legal for G. By partitioning into four cases based on the spins of v\ and 
V2, we can write the corresponding weighted partition function Z w (G) and its derivative as 



Z W (G)' = (A' + (wi + w 2 )A)z^z% 2 + (C + wxC)z^ + (D' + w 2 D)z^ + B', (8) 



G-S. 



□ 



Z W (G) = Az^zp + Cz^ + Dz 2 W2 + B; 



(7) 
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for polynomials A, B, C, D in the remaining variables z%. Notice that in the notation of Lemma 3.6, 
A = Z+({v\,v 2 }). Similarly, denoting the activity at the merged vertex by z, we have the fol- 
lowing expressions for G\ : 



Z W {G\) 



Az m+W2 + B; 

(A' + ( Wl + w 2 )A)z w i +W2 + B', 



(9) 
(10) 



with A and B as defined above. Now consider any fixing of the activities such that \zi\ > 1 for 
i > 2. Since G satisfies the weighted strict Gauss-Lucas property, we get by setting z\ = z 2 in 
eq. (8) that the (univariate) polynomial 

{A' + {wi + w 2 )A)z Wl+w * + (C + Wl C)z Wl + (D' + w 2 D)z w * + B' 

in z has no zeros satisfying \z\ > 1. Also, we know from Lemma 3.6 that A' + (w\ + w 2 )A / 0. 
Thus, we must have that the product of the zeros, B' / {A 1 + (w\ + w 2 )A), satisfies 



B' 



A 1 + [wi + w 2 )A 



< 1. 



However, using eq. (10), this implies that if \zi\ > 1 for i > 2, then Z W [G\)' can be zero only if 
\z\ < 1, and hence G\ satisfies the weighted strict Gauss-Lucas property. □ 

Finally, we give the proof of Lemma 3.4. 

Proof of Lemma 3.4. Note that any legal set of weights for Gi can be obtained by adding one to 
the weight w\ of v\ in a legal set of weights w of G, and then assigning u an arbitrary weight 
Wq > 1. With a slight abuse of notation, we denote these related weight assignments (one on G 
and the other on G\) by the same letter w. We now partition the terms in Z W (G) based on the 
spin of v\ to get 

Z w {G) = Az^+B- 
Z W (G)' = (A' + wiA)z™ 1 +B', 

where w\ is the weight of v\ in G. Here, A, B are polynomials in the remaining variables Zi, 
and A is of the form Z+({v\}) in the notation of Lemma 3.6. We again assume < f5 < 1 and 
\zi\ > 1 for i > 1. We now consider G\. Denoting the activity at u by z, we can write 



Z w {Gi) =A((3 + z 

Z w (Gx)' = (A' + w\A){fi + z W0 )z^ +1 + A(P + (w + l)z w °)z^ +1 + B'(l + p z w °) + w pBz 



+ B(l + /3z 



Now suppose that G\ does not satisfy the weighted strict Gauss-Lucas property, and hence \ z\ 
and \z\\ are both also at least 1, but Z W (G\)' = 0. It follows from Theorem 3.1 that we then also 
have Z W (G{) = 0. We now proceed to derive a contradiction to the above observations. For 
convenience, we denote z™ 1+1 by y in what follows. 

Using Lemma 3.6, we know that A ^ and that A' + w\ A ^ for our setting of activities. By 
Theorem 2.2 applied to Z W (G) and the weighted strict Gauss-Lucas property applied to Z W (G)', 
we get 



< 1, 



and 



B' 



A' + w\A 



< 1. 



Also, since Z W (G\) = 0, we must have 



B1 + j3z wo 
~A f3 + z w o 



(11) 



(12) 



11 



Notice that y is well defined since A ^ 0, \z\ > 1 and (3 < 1. Further, since (3 < 1, either one of 
\z\ > 1, or \B\ < \A\ would imply that \y\ < 1, which is a contradiction to our assumption that 
\zi\ > 1 (since y = z™ 1+1 ). Thus, we must have 



1, and 



1. (13) 



Now, substituting the value of y from eq. (12) into Z W {G\)' = 0, we get 

B 1 + (3z W0 



B'{1 + (3z w ») + Pw Bz W0 = ({A' + wiA) ((3 + z w °) + A{fi + (w + l)z w °)) 



A (3 + z w o 



Dividing through by (A' + w\A)(l + f3z w °), setting c = A/ {A' + w\A) and rearranging terms, 
we get 

B' B ( ( z W0 1 

1 + C + WqC — + — — 1 



A' + wiA A\ \l3 + z w o l + fiz w ° 

= \ { 1 + c + w oc (25ft (^w) - l) } , since \z\ = 1. (14) 

Notice that these divisions are well defined since A' + wiA ^ 0, and (3 < 1 and \z\ = 1 implies 
\hat{l+l3z W0 ) ^ Oaswell. Note also that c is of the form \/{wi+d) where K(c') = U(A' /A) > 
by Lemma 3.6 and our earlier observations about A: it therefore follows that K(c) > 0. However, 
we then calculate that for \z\ = 1, the factor inside the braces in (14) has real part (and hence 
absolute value) at least 1. Using \B\ / \ A\ = 1 from (13), we then see that the right hand side of 
(14) always has absolute value at least 1, which gives us the required contradiction to (11). This 
shows that G\ satisfies the weighted strict Gauss-Lucas property. □ 



4 Hardness of computing the mean magnetization 

In this section, we use our extended Lee-Yang theorem (Theorem 1.1) to prove Theorems 1.2 
and 1.3 via reductions from the problem of computing the partition function of the Ising model, 
which is known to be #P-hard even for bounded degree graphs [10, 17]. More specifically, we 
will use the following #P-hardness result. 

Theorem 4.1 ([10, Theorem 1],[17, Theorem 5.1]). Fix (3 satisfying < /3 < 1. The problem of 
computing the partition function Zi(G, (3, 1) of the Ising model on connected graphs of fixed maximum 
degree A > 3 is #P-hard. 

For simplicity, we prove here a version of Theorem 1.2 without the bounded degree con- 
straint. The extension to bounded degree graphs requires some more work and is proved in 
Appendix B. 

Proof of Theorem 1.2. We assume A > 1, since the case A < 1 is symmetrical. For given < f3 < 1, 
suppose that we have an algorithm A which, given a connected graph G, outputs the mean 
magnetization M(G,(3,\) in polynomial time. Let G be a graph of n vertices. Notice that as a 
rational function in z, M(G, (3, z) is a ratio of the two polynomials, VZj(G, (3, z) and Zi(G,f3,z), 
which are both of degree n. Further, since G is connected, these polynomials are co-prime by 
Theorem 1.1. Thus, if we could efficiently evaluate M(G, j3, z) at In + 2 distinct points z using 
algorithm A, we could uniquely determine the coefficients of Zj(G, f3, z) by Theorem 2.4 (since 
we know that the constant term in Zj(G, (3, z) is 1). We could then determine Zj(G, (3, 1) in 
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polynomial time. Theorem 4.1 would then imply that computing the mean magnetization for 
the given values of the parameters (3 and A is #P-hard. 

In order to evaluate M(G, j3, z) at 2n+2 distinct values, we consider the graph G(k) obtained 
by attaching k new neighbors to each vertex of v. We then have 

Zj(G(k), p, A) = (1 + pX) nk Zj{G, (3, X k ), and (15) 

k\{l- p 2 ) 



M(G«,/3,A) = ^ + 



1 + 



1 + /3AX/3 + A) 



M(G,p,X k ), (16) 



where X k = A (J$x) Notice that when p < 1, all the A fe are distinct, and further, M(G, P, A fc ) 

can be easily determined given M(G(k), p, A). Therefore, we can evaluate M{G{k), P, A) for 
< k < 2n + 1 using the algorithm A, and then using eqs. (15) and (16) we can determine 
M(G, P, Afc) in polynomial time. Since these evaluations are at distinct points, the reduction is 
complete. □ 

Proof of Theorem 1.3. For a given p as specified in the theorem, suppose that there is a poly- 
nomial time algorithm A which, given a graph G of maximum degree A, and a value of A in 
unary, outputs the susceptibility x(G, P, A). Notice that as a rational function in z, p, z) 
is a ratio of the two polynomials Z/(G, P, z) ■ V 2 Z I (G, P, z) - {VZ I (G, /?, A)) 2 and Z/(G, P, z) 2 , 
which are both of degree In. Further, since G is connected, these polynomials are co-prime by 
Theorem 1.1. To see this, notice that any common complex zero of these two polynomials must 
be a common zero of Zj(G, P, A) and VZj(G, P, A), which is prohibited by Theorem 1.1. 

To complete the reduction, we notice that we can choose 4n + 2 distinct values of A in the 
interval (0, 1] all of which have a unary representation length of at most 5n. Thus, using A, we 
can efficiently evaluate x{G, P, z) at 4n + 2 distinct values of z. By Theorem 2.4 we can then 
use these evaluations to uniquely determine the coefficients of Zj(G, P, z) 2 (since we already 
know that the constant coefficient is 1), and hence, Zi(G, p, 1), in polynomial time. Because of 
Theorem 4.1, this implies that the problem of computing the susceptibility at the given value 
of P is #P-hard. □ 



5 Hardness of computing the average dimer count 

In this section we prove Theorem 1.4 by reducing the #P-hard problem #Monotone-2SAT to the 
problem of computing the average dimer count. The reduction is similar in structure to Valiant's 
original proof for the #P-hardness of the problem of counting perfect matchings. However, 
since we will need to do rational interpolation, we need the zeros of the partition function to 
be simple, so by Theorem 2.3 we will need to ensure that the graph appearing as the output of 
the reduction always has a Hamiltonian path. The formal properties satisfied by our reduction 
are stated in the following theorem. 

Theorem 5.1. There exists a -polynomial time algorithm A which, when given as input a Monotone 
2-SAT formula (f>, outputs a weighted graph G with the following properties: 

1. The weights in G are drawn from the set {1, 2, 3}. 

2. Suppose </> has v variables and p clauses. Then, given the total weight W of perfect matchings in 
G, the number of satisfying assignments of 4> can be determined in polynomial time from W, p, 
and v. 

3. G contains a Hamiltonian path. 
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We observe here that Valiant's reduction from #3-SAT [36] can be easily modified so that 
it satisfies properties 1 and 2. However, it is property 3 that is crucial for our purposes, since 
it allows the use of Theorem 2.3. We first show how Theorem 5.1 can be used to immediately 
prove a slightly weaker version of Theorem 1.4, which shows hardness only on general graphs. 
The proof showing hardness for bounded degree graphs can be found in Appendix B. 

Proof of Theorem 1.4. Fix any A > 0, and suppose that there exists a polynomial time algo- 
rithm B which, given a connected graph H, with edge weights in the set {1,2,3} outputs 
D(H, (7 e ) eg £, A). In the following, we suppress the dependence on edge weights (■j e ) e& E f° r 
clarity of notation. Given a Monotone 2-SAT formula <fi, we can then produce the graph G = 
A{4>) in polynomial time. Let n be the number of vertices in G. Since G contains a Hamiltonian 
path, Theorem 2.3 implies that Z]\j(G, z) and VZm(G, z) have no common zeros. Thus, being 
able to use algorithm B to evaluate D(G, z) (and hence U(G, z)) at 2n + 2 different values of z 
would allow us to uniquely determine the coefficients of Zm(G, z) in polynomial time by ra- 
tional interpolation (Theorem 2.4), since we already know that the coefficient of z n is 1. This 
would allow us to obtain W (which is the constant term in Zm(G, z)), and hence, by property 2, 
also the number of satisfying assignments of 4>, in polynomial time. This would show that the 
problem of computing V(G, A) is #P-hard (since #Monotone-2SAT is #P-hard [37]). 

However, B only allows us to evaluate U (G, z) at z = A. In order to "simulate" other values 
of A, we consider the graph G(k) obtained by attaching k new vertices to each vertex of G with 
unit weight edges. We then have 

Z M (G(k),X) = X nk Z M (G,X k ); (17) 

U(G(k),X) = nk + ^^U(G,Xk), (18) 
X z + k 

where A& = A + k/X. Thus, by choosing 2n + 2 different values of k, none of which is equal to 
A 2 , we can determine U(G, z) at 2n + 2 different values of z by running B on G(k) and using 
eq. (18). This completes the proof. □ 

In the rest of this section, we proceed to prove Theorem 5.1 in a sequence of steps. For 
simplicity, we will describe our reduction in terms of cycle covers in a directed graph rather 
than perfect matchings in an undirected graph (this also allows us to directly compare our 
gadget construction with that of Valiant [36] at various steps). Given a weighted directed graph 
G = (V, E), we define the undirected bipartite graph Bip(G) = (V x {0, 1} , E') where the edge 
{(x, 0), (y, 1)} is in E' with weight 7 e if and only if (x, y) is an edge in E with the same weight. 
Note that a subset S C E forms a cycle cover of weight w in G if and only if the corresponding 
subset of edges S' = {{(x, 0), (y, 1)} \(x, y) G S} forms a perfect matching of weight win Bip(G). 
In particular, the total weight of all perfect matchings in Bip(G) is the same as the total weight 
of all cycle covers of G. 

Later, while arguing about the existence of Hamiltonian paths in graphs of the form Bip(G), 
we will find it convenient to use the following short-hand notation for simple paths in the graph 
Bip(G) in terms of the edges and vertices of G. Consider any simple path (x\, 1), (x2, 0), (23, 1), 
(a?4, 0), . . . (xi, 1), where we have assumed for simplicity that I is odd. The edges correspond- 
ing to this path in G are x\ <— xi, X2 — > x$, %s X4, . . . xi_\ — > x/. Notice that alternate 
edges are traversed in reverse in this representation. The path can therefore be represented 
as x\ <— X2 — > X3 <— X4 — > . . . — > xi. Similarly for a path starting on the other side, say 
(a?l, 0), (X2, 1), (xs, 0), (x4, 1), we have the representation x\ — > X2 <— £3 — > X4. Notice that a 
path P2 starting at a vertex v in this notation can be appended to a path p\ ending at v if and 
only if the arrows at v in p\ and P2 respectively are in opposite directions. We will refer to this 
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notation as the alternating path representation. Further, given an alternating path representa- 
tion of a path, we will refer to edges going right (such as x\ — > X2 in the last example) as forward 
edges, and edges going left (such as xi <— x% in the above example) as backward edges. 

5.1 Overview of the reduction 

We now look at the basic structure of our reduction, which is an elaboration of Valiant's reduc- 
tion [36] as modified by Papadimitriou [30] and presented in [1]. Recall that given a Monotone 
2-SAT formula <p, the reduction needs to produce in polynomial time a directed graph G such 
that the number of satisfying assignments of <fi can be easily determined from the total weight 
of cycle covers of G, and such that Bip(G) has a Hamiltonian path. Our first step is to introduce 
a shared variable in all the clauses of <fi: this shared variable will be useful later in showing the 
existence of a Hamiltonian path through the gadget. 

Observation 5.2. Let <p = Af=i °i ^ e a Monotone 2-SAT formula with /i clauses, v variables, and s 
satisfying assignments. Let rbea variable not appearing in <j) and consider the 3-SAT formula 

M 

^A(rVcj). 
i=i 

The number of satisfying assignments of '(f)' is s' := 2 U + s. 

Notice that each clause in <p' has exactly three variables, and that the number of satisfying 
assignments of <j) can be easily determined given the number of satisfying assignments of <j)' . 

We start the construction of G by creating a separate variable gadget (see Figure 1) for each of 
the variables r, x\, X2, ■ ■ ., x v occurring in <p' . This gadget has an external dotted edge for each 
appearance of the variable in the formula, and is designed so that any cycle cover must either 
use all the dotted edges in a particular gadget, or none of them. 

6* £ 

t -ft -ft' 

Figure 1: Variable gadget 

As done in Valiant's reduction, we then introduce a clause gadget (see Figure 2) for each 
clause in <p' . Each clause gadget has one external dotted edge for each literal in the clause, and 
is designed so that no cycle cover can include all the dotted edges; and so that for any other 
subset of the dotted edges, there is exactly one cycle cover including all the edges in the subset 
and no others. For each clause gadget, we label each of the three dotted edges in the gadget 
with one of the three literals appearing in the clause. However, in this step, we ensure that in 
each gadget the b — > c dotted edge is the one labeled with the literal r, since this is needed 
to show that the final construction has a Hamiltonian path. We now "pair" each dotted edge 
appearing in a clause gadget with a dotted edge corresponding to the same literal in a variable 
gadget, so that each dotted edge appears in exactly one pair. 

We first consider cycle covers which obey the constraint that they must choose exactly one 
edge from each such pair. We claim that the number of cycle covers satisfying this "pairing" 
constraint equals the number of satisfying assignments of </>'. To see this, we associate a truth 
assignment with every cycle cover by setting the variable v to true if the cycle cover uses all the 
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dotted edges in the variable gadget for v, and to false if it uses none of the dotted edges. Notice 
that because of the pairing constraint, a cycle cover is uniquely determined by specifying its 
assignment. Further, given the above properties of the clause gadget, exactly those cycle covers 
are permitted whose associated assignments are satisfying assignments of (/>'. 



We now enforce the "pairing" constraint referred to above using a gadget similar to Valiant's 
XOR-gadget. The XOR-gadget has two ports (labeled a and d), each of which admits one incom- 
ing and one outgoing edge (see Figure 3). To ensure the "pairing" constraint for a pair of dotted 
edges ei — > fi and e2 — > H, we replace them by the incoming-outgoing pair of a single XOR- 
gadget (see Figure 3b). The gadget has the property that after the replacement, the weight of 
every cycle cover which would have included exactly one of the two dotted edges e\ — > f\ and 
e2 — > /2 in the original graph gets multiplied by a factor of 2 (for each replacement made), while 
the weight of any cycle covers not satisfying the pairing constraint becomes (see Appendix D 
for a proof). The total weight of all cycle covers in the final graph so obtained is therefore 2 l s' ', 
where s' is the number of satisfying assignments and / is the total number of literals in <p' (since 
one XOR-gadget is needed to replace the pair of dotted edges for each literal). Further, replac- 
ing a pair of edges by a XOR gadget does not change the in-degree or out-degree of any vertex 
already present. 

Note that the XOR-gadget has edges of weight — 1, which are not permitted in the monomer- 
dimer model. This can be remedied by replacing the — 1 weight edges by a large chain of edges 
(of length, say, m 2 where m is the number of edges in the original graph) of weight 2, with 
individual vertices in the chain having self loops (of weight 1). The total weight of cycle covers 

2 

in the new graph modulo 2 m + 1 then gives the total weight of cycle covers in the original 
graph. 

This last step of replacing the — 1 edge by a long chain presents a challenge since we will 
need to include all the vertices in the chain in our Hamiltonian path (equivalently all — 1 weight 
edges must appear in the Hamiltonian path). For this reason, we cannot use Valiant's XOR- 
gadget directly. Our XOR-gadget, on the other hand, is such that the —1 weight edges can 
always be included in our Hamiltonian path. However, we have to be careful in the orientation 
of the XOR-gadgets in order to be able to construct a Hamiltonian path later: when replacing 
a pair of dotted edges one of which belongs to r's variable gadget, we orient the XOR-gadget 
so that the incoming edge at vertex a in the XOR-gadget comes from the variable gadget. At all 
other pairs, we orient the XOR-gadgets so that the incoming edge at the vertex a comes from a 
clause gadget. 




Figure 2: 3-SAT clause gadget 
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(a) XOR-gadget (b) Use of XOR-gadget 

Figure 3: Replacing dotted edges e\ — >■ fi and e2 — > H 

5.2 Analyzing the reduction 

We now proceed to analyze the output of the reduction to complete the proof of Theorem 5.1. 
The use of XOR-gadgets to enforce the "pairing" constraint as described above introduces a 
factor of 2 for each literal appearing in the clause, and therefore the total weight of cycle covers 
after this step is 2 3/ V < 2 6 ^. To get rid of the —1 weight edges in the XOR-gadgets, we replace 
each such edge by a chain of n = 6/x — 1 vertices with self-loops of weight 1 and connecting 
edges of weight 2. We call this final graph G. Since the initial total weight of cycle covers was at 
most 2 6/x , the weight of cycle covers in G (and hence the total-weight of all perfect matchings in 
Bip(G)), modulo 2 K+1 + 1, is exactly 2 3/ V. Since all steps in the construction of Bip(G) starting 
from 4> can be done in time polynomial in the representation size of <fi, this proves parts 1 and 
2 of Theorem 5.1. 

We now proceed to prove part 3, that is, that Bip(G) has a Hamiltonian path. We will use 
the alternating path notation described above in order to keep our discussion in terms of the 
vertices and edges of G, and we will call this representation of a Hamiltonian path in Bip(G) 
an alternating Hamiltonian path. In an alternating Hamiltonian path, each vertex of G is visited 
exactly twice, and the length of the alternating path between the two visits is odd. This is 
equivalent to saying that all vertices must appear exactly twice in an alternating Hamiltonian 
path, with all vertices except the first vertex in the path having one incoming forward edge, and 
one incoming backward edge. 

Our gadgets so far are designed to have alternating Hamiltonian paths which can be pieced 
together to form an alternating Hamiltonian path for G, and hence, we only need to list these 
paths and show how to stitch them together. We begin with alternating Hamiltonian paths in 
the clause gadget. 

Observation 5.3. The clause gadget in Figure 2 has the alternating Hamiltonian path 

c±- a <(— c — > b <r- a <(— b 
which uses all the dotted edges except theb c dotted edge. 

Recall that in the construction of the reduction, we ensured that the new variable r was 
associated with the b — > c dotted edge in each clause gadget. This will be used to connect the 
above alternating Hamiltonian path in different clause gadgets via connections to the variable 
gadget for r at the b — > c edge. Also, in the final construction, the dotted edges in the alternating 
Hamiltonian path will be replaced by detours into the associated XOR-gadget. 
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We now consider the XOR-gadget in Figure 3a. It turns out that in some cases, we will need 
to traverse the XOR-gadget partially, so that a path enters at a via a backward edge, uses the 
a — ^ d edge, and then leaves via a backward edge at d. In order to cover the rest of the vertices, 
we will then need to construct an alternating path that enters at a and leaves at d via forward 
edges, and covers all the vertices except a and d twice. Another complication with the XOR- 
gadget is the presence of two — 1 weight edges which need to be replaced with chains of vertices 
with self loops. However, this will not be a problem if we can ensure that both of the alternating 
paths described above use both the — 1 weight edges, since an edge in an alternating path can 
always be replaced by a chain of vertices with self-loops. We now show that, as we claimed 
above, our modified XOR-gadget satisfies all of these conditions. 

Observation 5.4. The XOR-gadget in Figure 3a has the alternating Hamiltonian path 

a <r- b — > b <r- a — > d -<r- c — > c ^ d. 

The gadget also has the following alternating path which enters at a and leaves from d using forward 
edges, but which does not otherwise visit these vertices: 

a <— 6—7-6^— c — >• c d. 

Moreover, both these paths use the — 1 weight edges b — > a and c — >■ c. 

Remark. Since the XOR-gadget is connected to variable gadgets (except those for variable r) 
in G via an outgoing edge at a and an incoming edge at d, it will be possible to replace the 
a — > d edge in the alternating Hamiltonian path above by a detour into the connected variable 
gadget when constructing an alternating Hamiltonian path in G. Similarly, it will be possible 
to use the a — > d edge as a replacement for the dotted edge in the variable gadget that was 
replaced by the XOR-gadget, at the cost of visiting the vertices a and d once. The role of the 
second alternating path is to visit the remaining vertices in a XOR-gadget which has already 
been partially traversed in this way. 

Remark. As pointed out above, it does not seem possible to include the two —1 weight edges 
in both the above alternating paths in Valiant's original construction. This necessitated the 
construction of our new XOR-gadget in which the — 1 weight edges are part of both the paths. 
The edges in our construction are the same as those in Valiant's construction, but the weights 
have been chosen differently. 

We now consider the variable gadget shown in Figure 1. We first work as if the dotted 
edges are present. In this case, for any vertex v in the gadget except the leftmost vertex, we 
can construct an alternating Hamiltonian path which covers all the vertices in the gadget, uses 
all the dotted edges except the one between v and its predecessor and can be appended to 
an alternating path that enters via a forward edge at v and leaves via a forward edge at v's 
predecessor. When the dotted edges are replaced by a XOR-gadget, this alternating path can 
still be traversed as described in the remarks following Observation 5.4, by instead following a 
forward edge to the d vertex of the XOR-gadget, following the a — > d edge in reverse, and then 
entering the variable gadget at the successor of v via the outgoing edge at the a vertex of the 
XOR-gadget. It is for this reason that we enforced above the condition that when a XOR-gadget 
is connected to a variable gadget for a variable other that r, it is oriented so that the incoming 
external edge at its d vertex comes from the variable gadget. 

We now start constructing the alternating Hamiltonian path in G starting at the left-most 
vertex in the variable gadget for r. If all the dotted edges were present, this gadget is just a 
chain of vertices, and hence there is an alternating Hamiltonian path that covers all its vertices. 
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However, each dotted edge has been replaced by an outgoing edge to the a vertex and an in- 
coming edge from the d vertex of a XOR-gadget. Thus, instead of following the dotted edges, 
our alternating path will take a detour into the corresponding XOR-gadget, and after traversing 
several other vertices, return via its d vertex to visit the other vertices in the variable gadget for 
t. Thus, we need to show that these detours into the XOR-gadgets can be used to make the al- 
ternating path go through all the other vertices in G twice while respecting the required parity 
constraints. 

We consider one such detour. We suppose that the XOR-gadget in question connects to a 
clause gadget C for the clause r V v\ V v<i- While following the alternating path for the XOR- 
gadget, we bypass the a — > d edge of the XOR-gadget and instead take a detour into the c vertex 
of C. We then start following the alternating path in Observation 5.3. If the dotted edges were 
present, we would be able to complete an alternating path covering all vertices in C and then 
return via a forward edge from the b vertex of C into the d vertex of the XOR-gadget. We could 
then complete the alternating Hamiltonian path in the XOR-gadget, and return via a forward 
edge into the variable gadget for r. However, since the dotted edges have been replaced by 
XOR-gadgets, we would need to take detours into the XOR-gadgets replacing them. Suppose 
we are trying to replace the dotted edge c — > a, corresponding to the literal v\. At this point 
there can be two cases: 

Case 1. Suppose that the vertices of the variable gadget for v\ have still not been covered by 
our growing alternating Hamiltonian path. Consider the XOR-gadget X replacing the c — > a 
dotted edge of C. We consider the alternating Hamiltonian path in Observation 5.4 starting at 
the a vertex of X. We follow this path until we need to use the a — > d edge. At this point, we 
take a detour into the variable gadget for v\ via a forward edge at vertex a of X. The vertex 
u we connect to in the variable gadget cannot be a leftmost vertex, since its predecessor u' 
is connected to vertex d of X via a u' — > d edge. As discussed above, we will therefore get 
an alternating Hamiltonian path for the vertex gadget which will leave the gadget through 
the u' — > d edge (though this will end up using the a — > d edges in all other XOR-gadgets 
corresponding to occurrences of v\). We can then complete the alternating Hamiltonian path 
for X, and this gives us an alternating path starting with a backward edge at vertex a of X, 
ending with a backward edge at vertex d of X, and covering all vertices in X and the variable 
gadget of v\, while also using up the a — > d edge in XOR-gadgets corresponding to all other 
occurrences of v±. We then use the forward edge from vertex c of C to vertex a of X and the 
forward edge from vertex d of X to vertex a of C to replace the dotted c — > a edge by the above 
alternating path. 

Case 2. Suppose that the variable gadget for the vertex v\ has already been covered by our 
growing Hamiltonian path. Then, as seen in Case 1, in the XOR-gadget X corresponding to 
the c — y a edge, the vertex a has already been visited using a backward edge, while d has al- 
ready been visited via a forward edge. Consider the second alternating path in Observation 5.4. 
Traversing this alternating path from a to d will satisfy the remaining covering requirements 
for all the vertices in X. Thus, as in Case 1 above, we can replace the c — > a dotted edge in the 
alternating Hamiltonian path for C by an edge from vertex c of C to vertex a of X and an edge 
from vertex d of X to vertex a of X. As before, this ensures that the vertices of this XOR-gadget 
are covered while traversing the alternating Hamiltonian path for C. 

Observe that since each clause gadget is connected to the variable gadget for r, and since all 
other variable and XOR-gadgets are connected to at least one of the clause gadgets, the above 
alternating path eventually covers all of the individual gadgets. This completes the proof for 
the existence of the alternating Hamiltonian path in G, and hence the proof of Theorem 5.1. □ 
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6 Future Work 



This work leaves open the complexity of computing several other average quantities; the most 
pertinent of which is (d), the average size of cuts under the Ising measure. The obvious ap- 
proach of attempting rational interpolation over ft via an analog of our Theorem 1.1 does not 
directly work, since Lee-Yang theorems do not hold in the same generality for the f3 parame- 
ter. A related problem is the complexity of computing the susceptibility x a * fixed values of 
A (in particular, A = 1), where, again, analyzing the partition function and its derivatives as 
polynomials in /3 may prove helpful. 

Extensions of our results to the antiferromagnetic Ising model and the hard-core (weighted 
independent sets) model also remain open: again, our current approach would need to be 
modified, since Lee-Yang theorems do not in general hold for antiferromagnetic systems. In 
particular, the best known analog of the Lee-Yang theorem for the hard-core model, due to 
Chudnovsky and Seymour [15], works only for claw-free graphs. Similarly, the complexity of 
computing averages in spin systems with more than two spin values (such as the Potts model 
or proper colorings) remains open. 

Finally, we mention potential connections with the large literature on stability preserving 
operators. As indicated in Section 1.3, this field is usually concerned with operators that pre- 
serve the region in which the zeros of a polynomial lie. Our extended Lee-Yang Theorem (Theo- 
rem 1.1) is an apparently rare example in which the operator actually makes this region strictly 
smaller. We conjecture that there may be more applications of this phenomenon. 

Acknowledgments. We thank Aris Anagnostopoulos for several fruitful discussions which led 
to the formulation of the problems considered in this paper. We also thank Marc Thurley for 
pointing out some useful references. 
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Appendix 



A Proof of Theorem 3.1 

As indicated earlier, Theorem 3.1 is an easy corollary of the following result of Newman [28] 
(restated in our notation). 

Theorem A.l ([28, Theorem 3.2 and eq. 3.5]). Consider the ferromagnetic Ising model on a graph 
G = (V, E) on n vertices, where we allow the edge potentials also to be variable, with the condition that 
the edge potential j3 uv on every edge uv satisfies < (3 UV < 1. Let (z v ) veV be a collection of complex 
vertex activities such that \z v \ > lfor all v G V. Then, 

71 

R(M(G,((3 uv ) uv&E ,(z) v&v ))>-. 

We now proceed with the proof of Theorem 3.1. 

Proof of Theorem 3.1. Let w be any weight assignment (not necessarily legal) of positive integral 
weights to the vertices of G. Consider the graph H obtained from G by appending to each 
vertex v of G a chain C v of w(v) — 1 vertices. Further, we let the edge potential be (3 on all 
edges of H which were present in G, and < 7 < 1 on all the edges which are either part of 
some C v , or connect a vertex v to its associated chain C v . Let (y v ) veV be a set of vertex activities 
on V. Henceforth, we will drop the subscript and refer to this set of activities as (y v )- With a 
slight abuse of notation, we also denote by (y v ) the collection of activities on H such that for 
any vertex x in H such that x G C v for some v G V, we have y x = y v . 

Now consider any collection of activities such that \ y v \ > 1 for all v G V. From Theorem A. 1, 
we get that for any 7 G (0, 1), 

77 

K(M(#,{/3, 7 },M))>~ (19) 

Since \y v \ > 1 for all v, the Lee-Yang theorem 3 implies that both Zi(H, {/3, 7} , (y v )) and Z W (G, (3, (y v )) 
are non-zero (when < 7 < 1). We can therefore take the limit 7 — > in eq. (19) to get 

\ < Jm K (M(H, {/3, 7} , (y v ))) = » (f^f^) = * ( M (G, 0, (y v ))) • (20) 

We now take a sequence ((yi))^ =1 of activity assignments such that \y*\ > 1 for a\\£ and v and 
such that lim y^ = z v . Since we assume in the hypothesis of the theorem that Z w ( G , j3 , ( z v ) v g y ) / 
0, we can take the limit i — > 00 in eq. (20) to get 

I < Hm K (m (G,P, (yi))) = »(M(G?)) , 

which completes the proof. □ 

3 Although we stated Theorem 2.2 only for uniform edge potentials, Asano's proof [2] in fact supports our current 
conclusion with variable edge potentials and a possibly disconnected graph. 
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B Hardness with a smaller blowup in degree 



Recall that, in our proofs of Theorems 1.2 and 1.4 in Sections 4 and 5, we realized different 
values of A required for the interpolation by attaching k extra vertices to each vertex of G. This 
necessarily entails a large increase in the degree of G. In this section, we give an alternative 
way of realizing different values of A which entails an increase in degree of exactly one, and 
which therefore allows us to complete the proofs of the stronger, degree-bounded version of 
Theorems 1.2 and 1.4. 

We denote by Pk a path of k vertices. Let p^ (respectively, p^) be the partition function 
Zj(Pk, f3, A) restricted to configurations in which the leftmost vertex of Pk is fixed to spin '+' 

(respectively, '— ')• We also set := %. Similarly, we denote by yk the partition function 

Zm (Pk, A), where we assume that all edges in Pk have weight one, and suppress the dependence 
on edge weights for clarity of notation. Notice that p± = 1 and p^ = r% = y\ = A. We further 
define yo = 1. The following recurrence relations show that, for fixed f3 and A, p\,Pk~, r& and 
yk can be computed in time polynomial in k: 

pt = HPPk-i+pti); (2i) 

Pk = PPk-l+Pk-V ( 22 ) 

. P + r k -i , oa . 
1 + Pr k -i 

Vk = Ayfc-i + yk-2- (24) 

Notice that p^ -,pZ, rt and y k are all functions of A. We note that values of their derivatives 
with respect to A can also be computed in time polynomial in k via the following recurrence 
relations: pj = yo = 0,pf = r\ = y\ = 1, and 

Pfc =P*A + AC8pfc_ 1 +pti); (25) 

Pk =PPk-l+Pk-V ( 26 ) 

. _ Pk Pk ~ Pk Pk . / 97 \ 
'k — , _ 2 > \ A/ > 

(Pk) 

y' k = Vk-i + Ay fc _i + ijk-2- (28) 

Here, we use the dot notation for the derivative with respect to A. Using a simple induction, 
one can also show that when /3 < 1, r& > for all k. 

Now consider a connected graph G. For k > 1, we define G(k) as the graph obtained by 
attaching to each vertex v of G a different instance of the path Pk, such that v is connected to 
the "leftmost" vertex of Pk via an edge. Notice that the maximum degree of G(k) is one more 
than the maximum degree of G. We first consider the Ising model on the graphs G(k). We have 

Zj(G(k),(3, A) = ( Pk - +1 ) n Zi(G, (3, A fc ), (29) 

where n is the number of vertices in G and A& = r^+i. Notice that when < j3 < 1, the 
sequence A& is strictly increasing and greater than 1 (respectively, strictly decreasing and less 
than 1) when A > 1 (respectively, when A < 1): this follows from the observation that the 
right hand side of the recurrence (23) is a strictly increasing function of r ( ! c _ 1 , and that r2 > r\ 
(respectively, r2 < n) when A > 1 (respectively, when A < 1). We also have 

M(G(k),/3, A) = + ^±iM(G, 0, A fc ). (30) 

Pk+i r k+i 

We now complete the proof of Theorem 1.2. 
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Proof of Theorem 1.2. As in the partial proof in Section 4, we assume A > 1 (since the case A < 1 
is symmetrical) and suppose that we have a polynomial time algorithm A which, given a con- 
nected graph G of maximum degree at most A > 4, outputs the mean magnetization M(G,/3,X) 
in polynomial time. 

Now consider any connected graph G of maximum degree at most A — 1 > 3. As shown in 
the partial proof in Section 4, Theorem 1.1 implies that if we can efficiently evaluate M(G, (3, z) 
at 2n + 2 distinct values of z using our hypothetical algorithm A, we can uniquely determine 
the coefficients of Zj(G, /3, z), and hence also Zj(G, /3, 1), in polynomial time. In view of Theo- 
rem 4.1, this would imply that the problem of computing the mean magnetization in graphs of 
maximum degree at most A for parameter values /3 and A is #P-hard. 

In order to evaluate M(G, f3, z) at 2n + 2 distinct values, we instead evaluate M (G(k), (3, A) 
for 1 < k < 2n + 2 using our hypothetical algorithm A. Notice that this can be done since 
the construction of the G(k) (as given in this section) implies that they have maximum degrees 
which are at most one larger than the maximum degree of G. Using eqs. (21) to (27) and (30), and 
the fact that > for all k, we can then determine M (G, f3, Xk) in polynomial time. Since Xk is 
a strictly increasing sequence, these evaluations are at distinct points, and hence the reduction 
is complete. □ 

We now consider the monomer-dimer model on the graphs G(k). We have 

Z M (G(k),\)=y%Z M (G,\ k ), (31) 

where n is the number of vertices in G and A& = yk+i/yk- We also have 

U(G(k), A) = n\t k + X(t k+1 - t k )U(G, X k ), (32) 

where t k = yk/Uk- It turns out that the sequence (\2k)k>o is strictly increasing and hence con- 
sists of distinct values, and further that t2fc+i — ^2k > for all k. This follows easily from the 
following explicit solutions for the yu and the Xk- 

£k+l _ ^k+1 £k+2 _ ^k+2 

where 

f = 1 ( A + x/A 2 + 4) > 0; v = l(x- \J\ 2 + A) < 0. 



2 V J '2 

Notice that t k +i — tk > for even k implies that for such k, we can determine U(G, Xk) given 
U(G(k), A), using equations (24), (28) and (32). We can now complete the proof of Theorem 1.4 
for the bounded degree case. 

Proof of Theorem 1 .4. As in the partial proof in Section 5, we fix any A > 0, and suppose that there 
exists a polynomial time algorithm B which, given a connected graph H with edge weights in 
the set {1, 2, 3}, and of maximum degree at most A > 5, outputs D(H, X) (recall that we are 
suppressing explicit dependence on the edge weights for clarity of notation). Given a Mono- 
tone 2-SAT formula cj>, we then produce the graph G = A{4>) in polynomial time. Notice that 
in the construction of G as given in Section 5.1, each vertex has degree at most 4 < A — 1: 
this corresponds to the maximum of the in-degrees and the out-degrees over all vertices in the 
directed version of the reduction. 

As argued in the partial proof in Section 5, Theorem 2.3 and the existence of a Hamilto- 
nian path in G together imply that if we could use algorithm B to evaluate D(G, z) (and hence 
U (G, z)) at 2n + 2 different values of z, then we can determine the number of satisfying assign- 
ments of 4> in polynomial time. This would in turn imply that computing D(G, A) for graphs of 
maximum degree at most A is #P-hard. 
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As before, in order to realize other values of A, we consider the graphs G(k) (as described 
in this section), for k = 0, 2, 4, . . . 4re + 4. Notice that the maximum degree of G(k) is one more 
than that of G, and hence is at most A. Further, as argued in the remarks following eqs. (31) 
and (32), these choices of k ensure that the values A& are distinct, and that U(G, A&) can be easily 
determined from U (G(k), A). We can therefore determine U (G, z) at In + 2 different values of 
z by running B on the G(k), as required. □ 



C Hardness for general two state feromagnetic spin systems and 
planar graphs 

We now show how to extend our results to general two-state ferromagnetic spin systems. Recall 
that a general two-state spin system [20] is parametrized by a (+, +) edge potential a\, a (— , — ) 
edge potential a 2 , and a vertex activity A. As before, given a graph G = (V, E), we define a 
probability distribution over the set of configurations a : V — > {+, — } via the weights ws{o) 
given by 

Ws {a) = \^a e 1 +{a) a e 2 {cT \ 

where e + (c) (respectively, e_(<r)) denotes the number of edges with '+' (respectively, '— ') spin 
on both end-points, while p(a) denotes the number of vertices with + spin. The partition func- 
tion Zs(G, «i, a 2 , A) and the magnetization M$(G, a\,a 2 , A) are given by 

Z s (G,ai,a 2 ,\) := ^ ws(<r); 

M,(g, a ,, n ,,A)^(p)- |7^»»w . 

Zs[G,ai,a 2 , A) 

Remark. The Ising model corresponds to the special case a% = a 2 = /3. 

It is well known that general two-state spin systems can be represented in terms of an Ising 
model in which the activity at each vertex depends upon the degree of the vertex [20]. In 
particular, if G is a A-regular graph then all vertex activities in the equivalent Ising model are 
the same, and one has 

ws(cr) = a 2 E){ wi{(j) (33) 

where the Ising model has an edge potential j3 = 1 / ^/aic^ and a vertex activity A' = A(ai/«2) A ^ 2 - 
A two-spin system is called ferromagnetic if the above translation produces a ferromagnetic Ising 
model, that is, when a\a 2 > 1. 

However, the above translation does not allow us to directly translate our hardness result 
for the ferromagnetic Ising model, since our results were not derived for regular graphs. We 
will instead do a reduction similar to the ones done in our earlier proof, but starting from the 
following somewhat stronger hardness result for the partition function. 

Remark. In this section, we allow graphs to have parallel edges and self-loops. In computing 
the degree of a vertex, each self-loop is counted twice (since it is incident twice on the vertex) 
and each parallel edge is counted separately. As observed in the remark following the proof 
of Theorem 3.2 in Section 3, our extended Lee-Yang theorem (Theorem 1.1) holds also in this 
setting. 

Theorem C.l ([13, Theorem 1]). Fix a\, a 2 > Qwitha\a 2 > land A > 3. The problem of computing 
the partition function Zg(G, ai,a 2 , 1) on A-regular graphs is #P-hard. 
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We will then prove the following theorem. We will also show later in this section that the 
theorem can be strengthened so that the #P-hardness holds even when the input is restricted to 
planar graphs. 

Theorem C.2. Fix a\, A > and A > 4 such that a\a2 > 1 . The problem of computing the 
magnetization Ms(G, «i, «2, A) on connected graphs of degree at most A z's #P-hard, except when a\ = 
a.2 and A = 1, in which case it can be solved in polynomial time. 

Remark. Notice that when 0102 = 1, the problem reduces to the case of a graph consisting of 
isolated vertices, and hence can be solved in polynomial time. Similarly, in the case a\ = a.2 
and A = 1, the two spins are symmetric, and the magnetization is therefore n/2, where n is the 
number of vertices in G. 

Before proceeding with the proof of Theorem C.2, we will need to analyze the model on 
graphs G(k) defined in Appendix B. As before, we begin by analyzing the model on the path 
P{k). We denote by p k (respectively, p k ) the partition function Z${Pk, «i, «i, A) restricted to 
configurations in which the leftmost vertex is fixed to be '+' (respectively, '— '). We also define 
the ratio r k =p\lp k - Similarly, we denote by m\ (respectively, m k ) the average magnetization 
of the path P k conditioned on the leftmost vertex being fixed to '+' (respectively, '— ')• We have 
pj = r\ = A,p^~ = 1 and = 1, = 0, and the following recurrences for k > 1: 

Pk =A(«iPfe-i+Pfc-l) (34) 

P k =a 2 p k _ l +pt_ l (35) 
air fc _i + 1 

r k = A ■ (36) 

a.2 + r fc _i 

mt = 1+ a ^t-iPU+rn^ 1 p^ 1 (37) 
aiP fe _i+P fc -i 

- a ^ m k-lPk-l + m k-lPk-l ,oox 

m k = _ (38) 

Under the condition a\a2 > 1, one can prove using a simple induction that for all k > 1, 
m k~ m k > 0/ and that when (ai — 1)A— («2 — 1) > (respectively, when (ai — 1) A— (02 — 1) < 0), 
the r k form a strictly increasing (respectively, strictly decreasing) sequence, and hence are all 
distinct. 

For reasons that will become clear shortly, we need the r k to be distinct, and hence we will 
need to handle the remaining case (a\ — 1)A — («2 — 1) = specially. We observe that unless 
a\ = a2 = 1, orai = 02 and A = 1, both of which are excluded in the hypothesis of the theorem, 
we cannot have both [a\ — 1)A — (02 — 1) = and {a\ — 1)A — {a\ — 1) = 0. To take advantage 
of this, we will modify P(k) by replacing each edge in P(k) by two parallel edges. We call 
the resulting graph P{k)' , and again define the quantities ptiPk- r k,mt anc ^ m k recursion 
on P{k)' . Notice that for k = 1, these quantities are the same as those for P(k); however for 
k > 2, we now need to modify the recurrences above by replacing a\ and «2 by a\ and a\ 
respectively. As before, we have — m k > for all A: > 1. Further, by our observation, the 
rfc form a strictly monotone sequence. Thus, in the case (ai — 1)A — (02 — 1) = 0, we redefine 
G(k) to use the paths P{k)' in place of P{k). In what follows, we will assume that G(k) are 
appropriately defined taking into account the values of A, a\ and 02, and will not explicitly 
keep track of the above modification. Notice that the maximum degree of G(k) is still at most 
max(A + 1, 3), where A is the maximum degree of G. 
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Given the above definition of G(k), we have the relations 



Z s (G(k),a 1 ,a 2 , X) = (a 2 p k + p^) n Z s {G, ai, a 2 , A fc ) 
M s (G(K),ai,a 2 , A) = nm k + (rri£ - m k ~)M s (G(k), ai,a 2 , X k ) 



(39) 
(40) 



where A& = X{a\rf t + l)/{a 2 + r k ). Since the r k form a strictly monotone sequence, it follows 
that (since a\a 2 > 1) so do the A&. In particular, all the Xk are distinct. 

Proof of Theorem C.2. Proceeding as in the proof of Theorem 4.1, we fix any A > and a\ , a 2 > 
satisfying a\a\ > 1, and suppose that there exists a polynomial time algorithm B which, given 
a connected graph H of maximum degree at most A > 4, outputs M$(G, a±, a 2 , A). 

Now consider any connected regular graph G = (V,E) of degree d := A — 1 > 3 on n 
vertices. From the translation in eq. (33), we see that for any A > 0, 



where (3 = \j ^Ja\a. 2 < 1. Theorem 2.4 along with our main Theorem 1.1 then implies that if 
we can efficiently evaluate M${G, a\ : a 2 , z) at 2n + 2 distinct values of z using our hypothetical 
algorithm B, we can uniquely determine the coefficients of Z${G, a\,a 2 , z), and hence also the 
value of Zs{G, ai, a 2 , 1), in polynomial time. In view of Theorem C.l, this would imply that 
the problem of computing the mean magnetization in graphs of maximum degree at most A 
for parameter values a\,a 2 and A is #P-hard. 

In order to evaluate Ms(G,a\,a 2 , z) at 2n + 2 distinct values of z, we instead compute 
Ms(G(k),ai,a 2 , A), for 1 < k < 2n + 2, using our hypothetical algorithm B. Notice that this 
can be done since the construction of the G(k) implies that they have maximum degrees which 
are at most one larger than the maximum degree of G. Using eqs. (34) to (38) and (40), and the 
fact that m k — m k > for all k, we can then determine M$(G, a\,a 2 , Xk) in polynomial time. 
Since A^ is a strictly monotone sequence as shown in the discussion above, these evaluations 
are at distinct values of z, and hence the reduction is complete. □ 

Extension to planar graphs 

Cai and Kowalczyk [13] also proved the following planar graph version of Theorem C.l. 

Theorem C.3 ([13, Theorem 1]). Fix ai, a 2 > with a\a 2 > 1, a / a 2 and A > 3. The problem of 
computing the partition function Zs(G, a\,a 2 , 1) on planar A-regular graphs is #P-hard. 

In order to extend Theorem C.l to planar A-regular graphs, we consider the cases qi ^ a 2 
and a\ = a 2 = a separately. In case a\ / a 2 , we proceed exactly as in the proof of Theorem C.2 
given above, except that we start with a planar ci-regular graph G in the reduction, and use 
Theorem C.3 instead of Theorem C.l as our starting hardness result. Since G is planar, so are 
the G{k), and hence we see that computing Ms {H, a\,a 2 , A) on planar graphs H, for a>i , a 2 and 
A satisfying the condition a\ ^ a 2 in addition to the conditions of Theorem C.2 is #P-hard. 

We now turn to the case a\ = a 2 = a > 1 (with A ^ 1). In this case, we start with 
the fact that computing Zs(G, 2a, ^, 1) on planar A-regular graphs is #-P hard (this is a direct 
corollary of Theorem C.3) . We again proceed exactly as in proof of Theorem C.2, starting with 
an arbitrary planar <i-regular graph G, and noting that the G(k) are planar too. Notice that the 
proof then shows that assuming the existence of a polynomial time algorithm to compute the 
magnetization in planar graphs of degree at most d + 1, we can evaluate the coefficients of the 
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polynomial Z S (G, a, a, z), and hence also the quantity Z$(G, a, a, 2 d ). However we then use 
the translation to the Ising model given above to see that 

Z S (G, a, a, 2 d ) = (|) Z s (g, 2a, |, l) , 

which shows that we can also evaluate Z$(G, 2a, ^, 1) in polynomial time. This establishes the 
#P-hardness in the remaining case ol\ = «2 (with A ^ 1). 

We thus see that in Theorem C.2, the input graphs can be restricted to be planar, and the 
same hardness result still holds. 

D XOR-gadget 

Claim D.l. The total weight of cycle covers of the XOR-gadget in Figure 3a is 2 when either 

• a is connected to an external incoming edge and d is connected to an external outgoing edge; 

or 

• a is connected to an external outgoing edge and d is connected to an external incoming edge. 

For all other external connections of a and d, the total weight of cycle covers of the gadget is 0. 

Proof. When the total number of incoming external edges at a and d is not equal to the total 
number of outgoing external edges, the XOR-gadget cannot admit a cycle cover due to parity 
constraints, and thus, the total weight of all cycle covers in these cases is trivially zero. A simple 
way to see this is that a cycle cover corresponds to a perfect matching in the natural undirected 
bipartite representation of the gadget discussed above. When the numbers of external incoming 
and outgoing edges are not equal, the bipartite graph does not remain balanced and hence 
cannot have a perfect matching. For all other configurations, in which the number of external 
incoming and outgoing edges are equal, the weights of all cycle covers can be shown to have 
the claimed value by exhaustive enumeration. □ 
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